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We discuss the theory of mixtures of Bosonic and Fermionic atoms in periodic potentials at zero temperature. 
We derive a general Bose-Fermi Hubbard Hamiltonian in a one-dimensional optical lattice with a superimposed 
harmonic trapping potential. We study the conditions for linear stability of the mixture and derive a mean field 
criterion for the onset of a Bosonic superfluid transition. We investigate the ground state properties of the 
mixture in the Gutzwiller formulation of mean field theory, and present numerical studies of finite systems. The 
Bosonic and Fermionic density distributions and the onset of quantum phase transitions to demixing and to a 
Bosonic Mott-insulator are studied as a function of the lattice potential strength. The existence is predicted of a 
disordered phase for mixtures loaded in very deep lattices. Such a disordered phase possessing many degenerate 
or quasi-degenerate ground states is related to a breaking of the mirror symmetry in the lattice. 
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I. INTRODUCTION 

Recent spectacular progress in the manipulation of neutral 
atoms in optical lattices (ll 0, |3|] has opened the way to the 
simulation of complex quantum systems of condensed mat- 
ter physics, such as high-Tc superconductors. Hall systems, 
and superfluid ^He, by means of atomic systems with per- 
fectly controllable physical parameters Optical lattices 
are stable periodic arrays of microscopic potentials created 
by the interference patterns of intersecting laser beams 01 ■ 
Atoms can be confined to different lattice sites, and by vary- 
ing the strength of the periodic potential it is possible to tune 
the interatomic interactions with great precision. They can 
be enhanced well into regimes of strong correlation, even in 
the dilute limit. The transition to a strong coupling regime 
can be realized by increasing the depth of the lattice poten- 
tial wells, a quantity that is directly proportional to the in- 
tensity of the laser light. This is an experimental parameter 
that can be controlled with great accuracy. For this reason, 
besides the fundamental interest for the investigation of quan- 
tumphase transitions 0] and other basic quantum phenomena 
HIESEIIIIEQ. optical lattices have become an impor- 
tant practical tool for applications, ranging from laser co oling 
fl4ll to quantum control and informationprocessing fTMfT^ . 
and quantum computation Ill7i.il8i.il9i.l20ll2lll . 

The theory of neutral Bosonic atoms in optical lattices has 
been developed f?] by assuming that the atoms are confined 
to the lowest Bloch band of the periodic potential. It can then 
be shown that the system is effectively described by a single- 
band Bose-Hubbard model Hamiltonian 153] . In such a model 
the superfluid-insulator transition is predicted to occur when 
the on site Boson-Boson interaction energy becomes com- 
parable to the hopping energy between adjacent lattice sites. 
This situation can be experimentally achieved by increasing 
the strength of the lattice potential, which results in a strong 
suppression of the kinetic (hopping) energy term. In this way, 
the superfluid-Mott-insulator quantum phase transition has 
been realized by loading an ultracold atomic Bose-Einstein 
condensate in an optical lattice 



The present paper is concerned with the study of dilute mix- 
tures of interacting Bosonic and Fermionic neutral atoms sub- 
ject to an optical lattice and a superimposed trapping harmonic 
potential at zero temperature. We assume the fermions to be 
identical (for instance spin-polarized in a magnetic trap), so 
that there are only s-wave Boson-Boson and Fermion-Boson 
contact interactions present. We construct an effective single- 
band Bose-Fermi Hubbard Hamiltonian, and we determine 
the ground state energy and on site density distributions in 
different mean field approximations of increasing complex- 
ity. Our main aim at this level of description is to determine 
the basic ground state properties of the mixture and to study 
how the Bosonic superfluid-insulator transition is influenced 
by the presence of the Fermions. Besides the study of the lat- 
ter issue, and the assessment of the properties of linear stabil- 
ity of the system against Boson-Fermion demixing, a remark- 
able finding of our analysis is that a quantum binary mixture 
loaded in a very deep optical lattice allows for a disordered 
phase of very many degenerate or quasi-degenerate ground 
states separated by very high potential energy barriers. In the 
limit of very large lattice potential strengths the basic mirror 
symmetry of the optical lattice is broken. 

The plan of the paper is as follows: In Section|II|we set the 
notation and derive the Bose-Fermi Hubbard model Hamilto- 
nian. We then discuss the range of validity of the approxima- 
tions and the assumptions used in the derivation. In Section 
mil we introduce some basic mean field descriptions to study 
the properties of stability of the mixture against phase separa- 
tion, and we provide a simple analytical criterion for the onset 
of a superfluid phase for the Bosons starting from a Mott- 
insulating ground state. 

In Section II VI we present numerical simulations for a 
small number of particles in the framework of the Gutzwiller 
variational ansatz. Two important cases have to be distin- 
guished: Boson-Fermion repulsive or attractive interaction 
(the Boson-Boson interaction is taken to be always repulsive). 
In the first instance, one can observe a continuous transition 
to a complete demixing of the Fermions from the Bosons, 
and to a Mott-insulating phase of the Bosons, as the strength 
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of the lattice potential is increased. In the case of attractive 
Boson-Fermion interactions there is no transition to demix- 
ing or to collapse of the mixture, while one still observes a 
Mott-insulating transition for the Bosons. By studying the 
behavior of the superfluid order parameter we show that in 
both cases the transition takes place at the same critical value 
of the lattice potential strength. Moreover, due to the strong 
attraction, the Fermions and the Bosons tend to form together 
ordered block-crystalline structures at the center of the trap. 

In Section|y]we present a numerical analysis that, although 
constrained to a small number of particles (five bosons and 
five fermions), seems to indicate the existence of a rich struc- 
ture of degenerate energy minima for large enough values of 
the lattice potential strength. Such an energy landscape sug- 
gests the possible existence of disordered phases of the mix- 
ture due to the delicate interplay between the different physi- 
cal parameters (Boson hopping, Fermion hopping. Boson on 
site energy, Boson-Fermion on site energy), the lattice depth, 
and the symmetries of the problem. Finally, a summary and 
an outlook to future research are shortly discussed in Section 

EH 



II. MODEL HAMILTONIAN 

We start by introducing the Hamiltonian for a Bose-Fermi 
mixture loaded into optical lattice potentials and confined by 
additional, slowly varying, external (harmonic) trapping po- 
tentials. It is given by 



H = Tb+Tf + Vb + Vf + Wbb + Wbf, (1) 

Tb = -Jdh^^{r)^^v), (2) 
Tf = -j£v¥{v)^i>{v), (3) 



where 



represent the Boson and Fermion kinetic energies, respec- 
tively, while 

Wbb - il^!!^!^ I d3r$t(i.)$t(r)$(r)$(r), (4) 



Wbf 



2 niB 
2TT^fPaBF 
mil 



d3j.<|t(i.)^t(i.)^(r)$(r) , (5) 



denote the Boson-Boson and the Fermion-Boson contact in- 
teraction energies. They are parametrized by the Boson- 
Boson and the Fermion-Boson s-wave scattering lengths 
ubb and gbf, respectively, and by the Boson mass tub and 
the reduced mass mji = mBruF / {ruB + mp), where niF 
denotes the Fermion mass. The potential energies 

Vb = j Sv¥{v){Vb{v)+Pb{v))^v), (6) 
Vf ^ j d^Y¥{r){VF{Y)+PF{Y))^{r), (7) 



are due to the trapping and lattice potentials. We consider pure 
magnetic trapping for Bosons and Fermions, so that Fermions 
are spin-polarized and their s-wave interaction energy Wf f 
can be neglected: Wff = 0. In the subsequent analysis 
we will consider the harmonic approximation of a typical 
quadrupolar magnetic field with strong anisotropy in the trans- 
verse directions y and z, i.e.. 



and 



VB{v)^mBLol{x^ + yy^ + yz^)/2, (8) 



VFiy) ~ ruFLol^x'^ + X^y^ + X^z'^)/2 , (9) 



where A ^ 1 is the anisotropy parameter. Moreover, if we 
assume trapping in the same magnetic state for the Bosons 
and the Fermions, then the trapping frequencies are related 
according to lof/lob = {mB/mFY/^, so that the two poten- 
tials coincide: Vb{y) = Vf{y)- The ground-state harmonic 
oscillator lengths, however, are different due to the different 
masses, and also differ for the a--direction on the one hand 
and the y and z-directions on the other hand: 



^11 

'■B/F 



= J^/imB/F^JB/F) 



(10) 



in the x direction, and H 



rections. We next consider a lattice structure for the Bosons 
and the Fermions in the x-direction, associated to the cor- 
responding Bosonic and Fermionic one-dimensional optical 
lattice potentials Pb{x) and Pf{x) are 



Pb{x) 

PFix) 



V^sm^{TTx/a) 
V^sm^{TTx/a) 



(11) 



where a is the lattice spacing associated to the wave vector 
k ~ 7r/a of the standing laser light. If the lattice poten- 
tials are produced by a far off-resonant laser for both species, 
the lattice potential strengths are equal for both Fermions and 
Bosons: Vp = Vg = Vq, and the two optical lattices coincide 
exactly. This is the situation we will always consider in the 
following. 

In the presence of a strong optical lattice and a sufficiently 
shallow external confinement in the x direction, the field op- 
erators can be expanded in terms of the single-particle Wan- 
nier functions localized at each lattice site .x,. Further, the 
typical interaction energies involved are normally not strong 
enough in order to excite higher vibrational states, and we can 
retain only the lowest vibrational state in each lattice potential 
well both for Bosons and Fermions (single-band approxima- 
tion). In case of stronger external confinements, or interac- 
tions, one should include higher Bloch bands as well in the 
expansion of the field operators, a case we do not consider 
in the present context. In the harmonic approximation, the 
Wannier functions w{r) factorize in the product of harmonic 
oscillator states in each direction, with the trapping potential 
almost constant between adjacent lattice sites. We then have 

<l(r) = Y.a,w^{x-x,)wf{y)wfiz), (12) 
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*(r) 



(13) 



where and bi are respectively the Bosonic and Fermionic 
annihilation operators at the i-th lattice site, Xi ~ ia, and the 
index i runs on positive and negative integers, the origin of the 
lattice being fixed at i = so that it coincides with the center 
(the minimum) of the external trapping potential. In each lat- 
tice potential well the Wannier local ground states for Bosons 
and Fermions are Gaussians in the harmonic approximation: 



,B/F 



(y) 



(^) 



exp 



B/F 



)l/2 



exp -zV2(^^/^)^ 



7,1/4 (£± 



B/F 



)l/2 



(14) 



(15) 



and 



exp -{x-x,)y2{£%/py 



where 



e%^p = a/[7r(Vo/E§/^y/^] 



(17) 



is the width of the harmonic oscillator potential wells at 
each lattice site, with E§ = {-kHY /2a^mB and Ep = 
{iiKf' /2a^mF being the Boson and Fermion recoil energies, 
respectively. 

In this paper we will consider the physical situation of 
very shallow trapping potentials, such that l\fp ^ aN^/F 
and consequently local density approximation (LDA) can be 
applied in the study of the ground-state properties of the 
system. Therefore, when exploiting the Wannier function 
expansions ( I12t and il3\ to map the full Hamiltonian Q 
into its lattice version, we discard all terms that are of or- 
der [aN B I p / t^^g I or of higher powers of it. Otherwise, 
nonlocal effects caused by the trapping potential, like site- 
dependent hopping terms, have to be considered f23ll. Finally, 
this approximation scheme leads to the following Hubbard- 
type Hamiltonian: 



H = -\Y. (-^sal+ia, + JFh\+^k + H. c 

■I 

+ "^T^f^l^ir^-D + UBFY. 



.(i).(i) 



F "-F 



(18) 



+ h {XiUB + Lo%/2) NB+h {XuJF + i^F/2) N, 



F ■ 



The first line in the above Bose-Fermi Hubbard Hamiltonian 
describes independent nearest-neighbor hopping of Bosons 
and Fermions, with amplitudes Jb and Jf, respectively. The 



terms in the second line describe Boson-Boson on site repul- 
sion (with Ub b > 0) and Boson-Fermion on site interaction. 
This interaction can be repulsive or attractive, depending on 
the sign of Ubf- The third line describes the energy offset at 
each lattice site due to the x component of the external trap- 
ping potentials Vg /^(r), and the last line contains the overall 
constant zero-point energy terms due to the y and z compo- 
nents of Vb/f{^) and to the lattice potential P{x). The on 
site interaction and offset energy terms are simple functions 
of the on site Boson and Fermion occupation number opera- 



tors 71^^ ~ a^cii and hp' — blhi, while the zero-point energy 
terms are proportional to the total particle number operators 
Nb = X^i ^l^i ^F = ^l^i- The frequency 

^l,F = m{li%/F?^B,F] (19) 

fixes the Bosonic and Fermionic harmonic oscillations in each 
lattice well. The relevant parameters entering in the Hamilto- 
nian are the on site values of the trapping harmonic potential 



(i) _ mB/F 2 2 
B/F ~ o ^B/F^i ' 



(20) 



the nearest-neighbor hopping amplitudes between adjacent 
sites Xi and Xi+i for Bosons and Fermions 



J, 



B/F 



dx w^^^ [x — Xi) 



2m B/F dx^ 



Vo siir I TT- 
a 



(21) 



the strength of the on site repulsion energy between two 
Bosonic atoms at the same lattice site 

inlfaBB [ i I Bi \\4. 

Ubb = / dx [w^ [x - Xi)) 



niB 



X / / dz{wf{z))\ (22) 



and the strength of the on site interaction energy (either repul- 
sive or attractive) between a Bosonic and a Fermionic atom at 
the same lattice site 



Ubf = 



2TTh?aBF 



dx [w^ (a; — Xi)wl^{x — Xi)] ^ (23) 



X J dy {y)]'j dz [u;f (z)u;f (z)]' . 

In typical situations we may neglect next-to-nearest neigh- 
bor hopping amplitudes and nearest-neighbor interaction cou- 
plings that are usually some orders of magnitude smaller, 
so that the Hamiltonian ( I18t provides a rather accurate 
model for the dynamics of a Bose-Fermi mixture with three- 
dimensional scattering in a one-dimensional periodic poten- 
tial. Terms involving nearest-neighbor interaction strengths 
and/or next-to-nearest neighbor hopping amplitudes can be- 
come relevant and need to be included, e.g., when consid- 
ering phonon exchange between Fermions, and this would 
lead to a Bose-Fermi analog of the so-called extended Hub- 
bard models. To evaluate estimates for the parameters enter- 
ing the Bose-Fermi Hubbard Hamiltonian ( I18> using Eqns. 
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(I14> . (I15> and (I16> . we will set the Boson recoil energy Eg = 
/i^TT^/ (2tobi^) as the unit of energy. We then introduce the 
dimensionless quantity Vb = Vq/E§, and, analogously, the 

dimensionless quantities tj bb^ U bf^ Vb 
We then have 



^'^ Vj,'\ Jb, and Jp- 



BB 



Ubf = 



V, 



8 ciBB a. -^1/4 



_ ■ , \ Ubf a 

3 i ^ mFj + 



(24) 



Fn'/' , (25) 



V, 



(l) _ TTlB 



(26) 



0.25 



0.15 




0.05 



Vo 



1 VJ) exp 



TT 



Jf 



- 1 Fo exp 



(27) 



(28) 



In FIG. [Owe show the dependencies of the these parameters 
on the potential strength Vq (compare also Ref. flai)- For ref- 
erence we have included as well the overlap integral 
{w{x — Xi)\w{x — Xi^i)) of adjacent Wannier functions. The 
overlap is negligible but for very small values of the poten- 
tial strength, confirming that terms of the order of the overlap 
integral can be neglected in the Hamiltonian. The Gaussian 
approximation holds rather well as can be seen by comparing 
the associated Bosonic hopping amplitude Jb with the one 
obtained by using the exact 1-D Mathieu equation 

Besides the conditions mentioned earlier, all the expres- 
sions derived in the present section are justified under the 
following circumstances. First of all, we must require that 
the two-body scattering processes are not influenced by the 
confinements, a condition that is guaranteed if the lengths of 
the confining and lattice potentials in all directions are much 
larger than the Boson-Boson and Fermion-Boson scattering 
lengths. Next, the single-band structure of the lattice Hamil- 
tonian is assured if the lattice spacing a is much greater than 
the harmonic confinements in each direction at all lattice sites. 
On the other hand, in this limit the harmonic approximation 
for the Wannier functions at each lattice well is automatically 
satisfied. Finally, as mentioned earlier, the assumption of a 
slowly varying confining potential such that LDA is applica- 
ble leads to the condition C-^b/f o-^b/f- We can sum- 
marize all the above conditions with the following chain of 
inequalities: 

{\aBF\,aBB 

(29) 

Our model is for some aspects unrealistic, since in present 
experimental situations the transverse confinements cannot be 
made very strong. Therefore a multi-band structure can ap- 
pear with several radial states being occupied, as reported in a 
recent experiment by the Florence group on Bose-Fermi mix- 
tures in a 1-D optical lattice ll24ll . 



FIG. 1: Top to bottom: the Fermion hopping amplitude Jf for 
mp /ruB = 0.5 in the Gaussian approximation; the Boson hopping 
amplitude Jb from the exact \-D Mathieu equation; the Boson hop- 
ping amplitude Jb in the Gaussian approximation; the Fermion hop- 
ping amplitude for mp/mB = 1-5 in the Gaussian approximation; 
and, for comparison, the overlap integral {w{x — Xi)\w{x — Xi+i)) 
of adjacent Wannier functions in the Gaussian approximation. All 
quantities are dimensionless. 



III. PHASE STABILITY AND THE SUPERFLUID 
TRANSITION 

In this section we investigate the zero temperature ground 
state properties of the system in a mean field approximation. 
In the following we will adopt a grand-canonical description 
through the Hamiltonian 



K = H - ijlbNb - ^ifNf , 



(30) 



where i^lb and ^p are the Bosonic and Fermionic chemical 
potentials. According to the Hohenberg-Kohn theorem, the 
ground state energy 



(*o|i^|*o) 



(31) 



is a functional of the on site Bosonic and Fermionic densities 



[a] hi 



and rip — {h\hi), where the expectation values 
are taken with respect to the ground state with state vector 
I'I'o). We decompose the functional E according to 



E — Eb ^" Ep -\- E 



BF 



- MS XI ^ ^^P^ 



(32) 



where Eb is the energy contribution depending only on the 

(i) 

Boson parameters Jb, Ubb, Vb , Ep is the energy depend- 
ing only on the Fermion parameters; and Ebf is the term 
due to Boson-Fermion interactions. We treat this latter term 
in mean field approximation: neglecting exchange correlation 
effects: 



Ebf = Ubf ^ n^BrSp' 



(33) 
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Exchange correlation effects have been recently studied for 
the case of homogeneous mixtures in the continuum ll25ll2^ . 
For the Fermion energy Ep^ we take the energy of the non- 
interacting homogeneous system and exploit local density ap- 
proximation (LDA) on it. 



2Jp 



^sii 



(34) 



This approximate description of the Fermions is well justified 
in the presence of a slowly varying trapping potential (so that 
LDA can be applied), when there are no direct interactions 
among the Fermions (as in our case), and moreover when one 
can neglect induced phonon-mediated self-interactions due 
to the presence of the Bosons. Therefore, in this situation, 
the nontrivial features of different quantum phases will regard 
only the Bosonic sector and not the Fermionic one. However, 
the presence of the Fermions will indirectly contribute to the 
properties of the different Bosonic phases, and this is the sub- 
ject that we will study in the following. 

In order to find an expression for the Boson energy Eb we 
will proceed in steps of increasing accuracy. First we per- 
form a very simple mean field analysis in two extreme lim- 
its: a completely superfluid Boson ground state and a totally 
Mott-insulating Boson ground state. In the latter case we will 
provide a simple criterion for stability of the mixture against 
demixing. Next, we will perform a perturbation expansion 
around the Mott-insulating Boson ground state to recover per- 
turbatively the phase boundary against transition to superflu- 
idity. Finally, in the next section, we will study the ground 
state properties of the mixture using a Gutzwiller ansatz for 
the Bosons capable of describing the intermediate regimes be- 
tween the insulating and superfluid Bosonic phases. 

We first consider the Bosons to be superfluid. In this regime 
the chemical potential and the number of particles in a homo- 
geneous system are related, to lowest order in Ubb, via 



Mb = Ubbit-o - 2Jb 



(35) 



where no is the density of condensed Bosons. Additionally, 
for very weak interaction uq « ub- Exploiting this result in 
LDA and using the mean field expression for the Bose-Fermi 
interaction energy we can then write for the inhomogeneous 
Bose-Fermi mixture at a given lattice site: 



Ubbu'-;^ = ^iB + 2Jb - V^^ - UBpn'-p 



(36) 



Next, we consider the case of a Mott-insulating Bosonic 
phase. To lowest order in Jb we neglect the kinetic term al- 
together. Then it is easily shown that the relation between 
the Bosonic chemical potential and the Bosonic density for a 
homogeneous system is given by 



fJ-B = UBBnB - UbbI'2 



(37) 



Exploiting LDA as before, we have in the inhomogeneous 
case at a given lattice site: 



Ubbu^^ ^^1B + Ubb/2 - 



Ubf-ti^f 



(38) 



Comparing Eqns. ( I36> and ( I38> . we observe the same behav- 
ior of the on site density profiles but for a constant correc- 
tion to the Boson chemical potential depending whether the 
Bosons are in a superfluid or in a Mott-insulating state. Fi- 
nally, differentiating the energy functional with respect to the 
on site populations of the Fermions, we determine the associ- 
ated density field and the set of coupled equations describing 
the ground state of the mixture at any lattice site. 



Ubbti-b 



^^'B 



V, 



Ubftl^f , 



(39) 



(40) 



where /i'^ is the proper expression of the Boson chemical 
potential according to whether the Bosons are in the Mott- 
insulating or superfluid regime. These equations are valid at 



a given lattice site i if ^'^ 



V, 



UBpny > 0, otherwise 



one must set n . 



0. On the other hand, if 



i^iF-V^'^ -UBFn^B^)/2JF <0 



(41) 



we must impose n^"* = at the given lattice site, while rip^ = 

1 must be imposed when {fiF ~ Vp^ — UBpn^B) I'^'^p ^■ 
These expressions are the lattice analogs of the Thomas- 
Fermi description of Boson-Fermion mixtures in the contin- 
uum. We remark that in the Mott-insulating regime the Bo- 
son on site populations rig must be rounded off to the integer 
closest to the solutions of Eqns. ( I39> -( I40> . 

In the Mott-insulating regime we can determine a criterion 
of linear stability against phase separation of the two species 
if we expand the energy functional E to second order in the 
small density variations Sn'^jp around the minimum provided 
by the solution of Eqns. j39>-j4m: 



5'E 




Ubb 



Ubf 



Ubf 2TrJFSm{Trn^p) 



(42) 



This quadratic form is positive at a given site i if and only if 



27rJFsin(7m^^)[/BS > U^f 



(43) 



and 27rJ_F sin(7rny ) + Ubb > 0. This last condition is al- 
ways satisfied for Ubb > and identical Fermions. If this is 
not the case for every site i, then the ground state is not stable 
against demixing. This result is similar to that recently ob- 
tained for a mixture of two different Boson species on a lattice 
l27il . which states that the mixture is stable if U1U2 > U12, 
where Ui and U2 are the Boson-Boson interaction strengths 
of species 1 and 2 respectively, and U12 is the interspecies 
coupling. The form of expression ( I43l l then suggests that the 
Pauli on site energy 27r Jj? sin(7rn^'') has the meaning of a 
density-dependent interaction strength. A similar correspon- 
dence was previously pointed out for homogeneous Bose- 
Fermi mixtures in the continuum 1291 . 
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Introducing a perturbation expansion with respect to Jb 
around the Mott-insulating ground state we can recover the 
zero-temperature phase transition to the superfluid phase. The 
reverse, i.e., to build a perturbative expansion in powers of 
Ubb around the superfluid ground state fails to describe the 
transition to a Mott insulator, as pointed out in Ref. 01 for the 
pure Bose case. We follow the procedure adopted in Ref. l27il 
for the two-component Boson mixture, with the due modifi- 
cations for the present case of a Boson-Fermion mixture, by 
treating the Bosonic kinetic (hopping) term as the perturbation 
with respect to the Bosonic Mott-insulating ground state. This 
scheme was first introduced for one-component Bose systems 
in Ref s. tSt, .29^ .3Q1 . We proceed by expanding the ground 
state energy with respect to the (local) Bosonic superfluid pa- 
rameter ■0'^') = [ai). At the phase boundary between a Mott 
insulator (MI) and a superfluid (SF) the expansion coefficients 
must vanish, yielding the following criterion for the onset of 
the transition to the (local) SF state: 



i7BB(27i^'-l)-2J£ 



uIb - ^uIb 



(2n« 



1) 



4J| 



1/2 



ii) 



Ubfi4^ 
- 1) - 2 Jb 



u'bb - 



(2n« 



1) + 4J| 



1/2 



(44) 



The minimum value of Ubb I Jb, where a MI phase can exist, 
is given by the condition 



C/bb/ Jb = 4n^) + 2 + 2V(2n 



1)2 - 1 



(45) 



conserving Gutzwiller ansatz in the determination of the en- 
ergy functional, while keeping the same approximations pre- 
viously introduced for the Boson-Fermion interaction and the 
Fermion energy, the total ground state energy reads 



E = Eb+Ep + Ubf^ 



^B ' 



(47) 



where the subsidiary conditions ensuring particle number con- 
servation are 



E 



{a\a^) ^ Nb, 



(48) 
(49) 



The Boson energy contribution is now 

Eb = -ijB(^^^(^+^)V 



C. c. 



(50) 



and the Bosonic observables are related to the Gutzwiller am- 
plitudes by 



= E (51) 

oo 

= {a\aia\a,)=Y,n'{fl:^)\ (52) 

oo 



and it involves the Fermionic sector indirectly through the de- 
pendence of ny on the Fermionic parameters and density 
distributions provided by Eqns. (I39ll-( l40t . Apart from this 
important modification, the phase diagram, at this level of ap- 
proximation, is analogous to that of a one-component Bose 
system. 



IV. NUMBER-CONSERVING GUTZWILLER ANSATZ 
AND NUMERICAL ANALYSIS 

The simplest ansatz for the Boson ground state being ca- 
pable of describing both the SF and the MI phases is the 
Gutzwiller Ansatz, which contains the mean field approxima- 
tions previously discussed as special cases. It consists of fac- 
torizing the amplitudes of superpositions of all possible Fock 
states consistent with a fixed number of Bosons Nb, in the 
following way L31i1 ; 



-|o>. 



(46) 



In this way correlations between the Bosons under particle 
exchange are included in a natural way. Using the number- 



Moreover, we must impose the natural constraints that 



n=0 

< < 1, 



(54) 
(55) 



for each lattice site i, reflecting the fact that the Gutzwiller 
amplitudes form a probability distribution for each lattice site, 
and that the on site Fermion occupation number cannot exceed 
one. 

To identify the ground state amounts to solving a con- 
strained optimization problem: one has to minimize the en- 
ergy functional i50\ subject to the constraints given by Eqns. 
( B^ and together with Eqns. and (|53}. We have 
solved the problem numerically for a small system of ten par- 
ticles (five Bosons and five Fermions). The first observation 
is that the optimization problem is not a convex optimization 
problem. Hence, one has to expect several local, "poorer" ex- 
trema in addition to the (not necessarily unique) global one. 
The numerical solution of this optimization problem has been 
performed first using a simulated annealing method l32il with 
an appropriate logarithmic annealing schedule. The quadratic 
constraints (|53 and ( I55> have been incorporated in a dynam- 
ical penalty formulation (see, e.g., Ref. 13311 ). Finally, for the 



FIG. 2: On site Bosonic densities for a Bose-Feimi repulsion ubf = 
0.04, as a function of the lattice potential strength. In this figure - as 
well as in the following figures - Vq runs from 1 to 8. 



local refinement the Nelder-Mead downhill simplex method 
ll3"4i1 has been applied. 

In FIG.|2]we show the change of the on site Bosonic den- 
sities with increasing lattice potential strength Vq for a sys- 
tem of five Bosons and five Fermions with moderate repulsive 
Boson-Fermion interaction. We note from FIG.IHthat as the 
strength of the lattice potential increases the Bosons go in a 
complete Mott-insulating phase, forming a block crystalline 
configuration around the center of the trap (which coincides 
with the origin of the optical lattice) with exactly one Boson 
per lattice site. The corresponding on site Fermionic densities 
are plotted in FIG. |3l From both figures we can see that, if 
Ubf > 0, by increasing the lattice potential strength the sys- 
tem eventually undergoes simultaneously a Boson MI transi- 
tion and complete phase separation, in accordance with Eqn. 
( l43l along with Eqns. (l24l - (l28t. 




FIG. 3: On site Fermionic densities for a Bose-Fermi repulsion 
asF = 0.04, as a function of the lattice potential strength. 

The local Bosonic superfluid parameter i/i'^'^ for the same 
physical situation is shown in FIG. |4] Although we are deal- 
ing with a finite system, we can already see a rather clear sig- 
nature of the onset of a phase transition to a Mott-insulator 
regime when the superfluid parameter suddenly drops to very 
low values at an approximate critical lattice potential strength 
Vo' ^ 7. 

We next consider the ground-state properties in the case 
of an attractive Boson-Fermion interaction. Because of the 



FIG. 4: The Bosonic superfluid on site order parameter for a Bose- 
Fermi repulsion asF = 0.04, as a function of the lattice potential 
strength. 




FIG. 5 : On site Bosonic densities for a Bose-Fermi attraction asF = 
— 0.04, as a function of the lattice potential strength. 



strong attraction with growing lattice depth, the Fermions fol- 
low the Bosons in building a sharp crystalline block around 
the center of the trap, as can be seen from FIGs. |5] and 
|6| We cannot expect in this case to observe a simultaneous 
mean field collapse like the one predicted for a trapped Bose- 
Fermi mixture in the continuum llH Is^l (for the effects be- 
yond mean field see IjtIi '). as this possibility is forbidden in a 
single-band approximation. Finally, we consider the behav- 
ior of the Bosonic superfluid on site parameter in the case of a 
Boson-Fermion attractive interaction. Comparing FIG.^lwith 




FIG. 6: On site Fermionic densities for a Bose-Fermi attraction 
Ubf = —0.04, as a function of the lattice potential strength. 
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FIG.0] we see that the transition to a Mott insulating phase for 
the Bosons takes place at the same lattice potential strength, ir- 
respectively of the repulsive or attractive nature of the Boson- 
Fermion interaction. This finding confirms the results of the 
mean field analysis presented in the previous Section. 




FIG. 7: The Bosonic superfluid on site order parameter for a Bose- 
Fermi attraction ubf = —0.04, as a function of the lattice potential 
strength. 



V. MIRROR SYMMETRY BREAKING AND TRANSITION 
TO DEGENERACY 

The above optimization problem associated with the con- 
strained minimization of the energy is not convex, hence there 
can be many local minima in addition to the global one. How- 
ever, even the ground state may be approximately or exactly 
degenerate. In fact, this is what happens in the case of Boson- 
Boson and Boson-Fermion repulsion for large values of the 
lattice potential strength Vq. As Vq grows, it becomes eventu- 
ally energetically more favorable for the bosons to be arranged 
in single-particle occupancy of the available sites around the 
center of the external trap. The Bosonic and Fermionic on- 
site occupation numbers can only assume the values or 1, 
and a definite Boson-Fermion symmetry is established in the 
Bose-Fermi Hubbard Hamiltonian assuming that the on site 
Fermionic and Bosonic trapping potentials coincide. 

A similar transmutation of Bosons into Fermions in strong 
optical lattices has been pointed out by Paredes and Cirac in 
a recent paper fllll . They consider a model of pure Bosons 
in an optical lattice and show that in the limit of very strong 
Boson-Boson on site interaction, the Bosonic operators can 
be mapped into Fermionic operators by means of the well 
known Jordan- Wigner transformation. Let us consider what 
happens in the case of a Boson-Fermion mixture. As the lat- 
tice strength grows, configurations of lowest energy that are 
mirror-symmetric with respect to the center of the lattice, like 
e.g. those of FIGs. 2 and 3, become approximately ener- 
getically equivalent to other symmetric configurations (e.g., 
a checkerboard of alternating Bosons and Fermions with one 
particle per lattice site), as well as to nonsymmetric configu- 
rations (like a succession of four Fermions followed by five 
Bosons and then a last Fermion, again with one particle per 
lattice site), and mirror symmetry breaking takes place. We 



may thus consider sequences of energy functionals with in- 
creasing lattice potential strengths Vq- For each value of Vq, 
one may identify a ground state. Then, the difference in en- 
ergy of this ground state to those states that can be obtained 
by interchanging the role of Fermions and Bosons will con- 
verge to zero as Vq grows. The Boson hopping contribution 
will become negligible, whereas the behavior of Ub b will en- 
force the mean Bosonic on site occupation number to be at 
most one. Hence, for each lattice site, the constraints on the 
Boson and Fermion occupation numbers become identical (at 
most one Boson or one Fermion per lattice site). Notice that 
the suppression of the hopping terms is exponential. More- 
over, since Vg'' = Vp'' for all lattice sites i, the larger the 
value of Vo, the more symmetric is the role of Bosons and 
Fermions. There are then many ground states that are de- 
generate in energy with respect to any permutation of lattice 
sites - as long as all particles are located around the mini- 
mum of the confining external potential V^^"* — Vp^'' = 0. 
These degenerate configurations will be given by all possi- 
ble symmetric and nonsymmetric Fermion and Boson distri- 
butions in a region around the center of the lattice, with ev- 
ery site of the region occupied by one and only one particle. 
Such possible configurations are for example checkerboard al- 
ternating patterns of Bosons and Fermions, or Mott Bosonic 
central configurations with Fermionic wings on the sides, or 
consecutive block crystalline arrangements of variable length 
of Bosons and Fermions. In brief, while the Hamiltonian for- 
mally retains its mirror symmetry under reflection of the lat- 
tice around its center, the degenerate ground states need not, 
and spontaneous mirror symmetry breaking occurs. At the 
same time complete Boson-Fermion exchange symmetry sets 
on. No ground state is a priori favored compared to any other: 
any random pattern of consecutive Bosons and Fermions lo- 
cated around the minimum of the external trapping potential 
is a legitimate ground state. FIG. |8] shows representative on 
site Bosonic densities in the regime of large values of Vq 
around Vq = 50 for the case of Boson-Boson and Boson- 
Fermion repulsion in a system composed of five Bosons and 
five Fermions: at each value of the lattice potential strength, 
a particular state is selected from the set of those with same 
energy. Each vanishing value of the on site Bosonic density 
means that exactly one Fermion has filled that particular lat- 
tice site. The large value chosen for Vq allows to clearly stress 
the random nature of the configuration patterns even for very 
small changes of the lattice potential strength, whereas degen- 
eracy and disorder can set in already at lower values of the lat- 
tice depth, depending on the tuning of the harmonic oscillator 
and scattering lengths (see below). The degenerate states are 
separated by energy barriers. The system is non-ergodic, and 
hysteresis should be observed: what particular state is chosen, 
depends on the exact mechanism of preparation of the system 
and of loading of the mixture into the optical lattice. 

The criterion for the onset of degeneracy and nonperiodic 
ground states in the bulk region around the center of the lattice 
and of the trapping potential is easily identified, by looking 
at the relative importance of the trapping on site energy with 
respect to the on site Boson or Fermion interaction energy. For 
instance, to allow for the Fermionic behavior of the Boson on 
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FIG. 8; The disordered pattern of Bosonic ground-state distribu- 
tions for repulsive Boson-Boson and Boson-Fermion interactions 
for large values of Vo around Vo = 50. 

site occupation numbers (either or 1) one must require that 
the energy is lower having one Boson at the edge of the bulk 
central region rather than having it sitting on top of another 
Boson at the center of the lattice: 

Ubb >Vb{i^ (Nb + Np)/2) . (56) 

The analogous condition for the Bose-Fermi on site interac- 
tion is; 

Ubf > Vb/f = {Nb + Nf)/2) . (57) 

For smaller values of Vq, the Boson hopping contribu- 
tion will become more and more important. A representa- 
tive situation of this intermediate regime is depicted in FIGs. 
12] and |3] here, the repulsion between Bosons and Fermions 
is strong enough to allow for phase separation, while the 
non-negligible hopping terms still favor configurations where 
Bosons have Bosons as nearest neighbors. The transition to 
degeneracy and disorder, exact in the limit of infinite lattice 
depth, is a novel and peculiar feature of Bose-Fermi mixtures 
and it should hold in general for any multicomponent Bose 
and/or Fermi dilute atomic system loaded in a deep optical lat- 
tice at zero temperature, provided that intercomponent inter- 
actions are repulsive and the on site confining potentials coin- 
cide for the different components. It clearly cannot take place 
in a single-component system, say a pure single-component 
Bose gas, where only a SF-MI transition occurs 0]. The 
rather complex and rich interplay between ordered and dis- 
ordered configurations of Bose-Fermi mixtures in very deep 
optical lattices will be considered in more detail elsewhere. 

VI. SUMMARY AND OUTLOOK 

In conclusion, we have studied the zero-temperature prop- 
erties of a mixture of weakly interacting gases of neutral 
Bosonic and Fermionic atoms loaded in one-dimensional op- 
tical lattices and confined by harmonic trapping potentials. 
We have derived a single-band Bose-Fermi Hubbard Hamil- 
tonian, and performed some mean field studies of the zero- 
temperature phase diagram. We have considered the case of a 



quasi-free Fermion sea acting on the Bosons, which have been 
treated in their full dynamical range. We have always worked 
in the approximation of s-wave Boson-Boson and Boson- 
Fermion contact interactions. According to the possible dif- 
ferent combinations of intraspecies and interspecies attractive 
and repulsive interactions, the system displays a rich phase 
structure, including the onset of a SF-MI transition in the Bo- 
son sector, and a simultaneous transition to demixing in the 
Boson-Fermion sector. The optical lattice potential plays a 
crucial role, allowing to tune the system into regimes of strong 
Boson-Boson and Boson-Fermion couplings as the lattice 
depth is increased. For very deep lattices the system displays a 
remarkable transition to a multiply degenerate phase in which 
all possible permutations of configurations with one Bosonic 
or Fermionic atom per site are legitimate ground states. The 
transition is related with breaking of the lattice mirror sym- 
metry for very large values of the lattice depth. This peculiar 
disordered pattern of degenerate ground-state configurations 
separated by very large barriers is somehow reminiscent of 
the behavior of classical disordered systems like glasses and 
spin glasses, but it takes place in a quantum system at zero 
temperature. 

The setting that has been investigated in detail in the present 
paper can be extended in various ways. Certainly a larger 
number of Bosons and Fermions have to be considered in or- 
der to obtain a more realistic description of the system. While 
our previous analytical findings are applicable to any num- 
bers of atoms, in order to extend the numerical calculations to 
larger numbers more powerful numerical methods have to be 
introduced. So far, Monte Carlo simulations with a fairly large 
number of particles have been carried out only for an inhomo- 
geneous Bose-Hubbard model jsH- The authors have also 
speculated that the qualitative phase diagram does not depend 
on the dimensionaUty of the system. 

In order to extend the work presented in this paper to 
the case of interacting Fermions, one may either allow for 
different Fermionic species in magnetic traps or for spin- 
unpolarized identical Fermions in optical traps. Exchange- 
coiTelation effects, that are already included in the Gutzwiller 
Ansatz for the Bosons, will then become important for the 
Fermions as well Isoll . In this way the Bose-Fermi interac- 
tion has to be consistently incorporated beyond the mean-field 
level. In the present paper the mean-field treatment of the 
Bose-Fermi interaction is consistent, since we only included 
the Fermions in a mean-field and LDA prescription. 

At the opposite extreme, one can consider the Fermions to 
be static impurities for the Bosons, as very recently discussed 
by Vignolo and collaborators 14^ . Finally, after the first ver- 
sion of the present paper had been submitted, a preprint ap- 
peared on the induced Boson-Boson interaction due to the 
Fermions for mixtures loaded in a 2-D lattice j4lll . 

Besides these fundamental theoretical aspects related to the 
theory of quantum phase transitions and the statistical me- 
chanics of complex systems, ultracold Bose-Fermi mixtures 
in an optical lattice qualify for potential applications in the 
physics of quantum information. As with systems involv- 
ing either Bosons or Fermions that have been studied so far 
mixtures could be used for the preparation 
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of multi-particle entangled states 1 1611 such as cluster states or 
certain instances of graph states 1431 . as well as for the imple- 
mentation of quantum gates. With Bosons and Fermions serv- 
ing two different purposes, Bose-Fermi mixtures could in fact 
allow for new possibilities of quantum information processing 
in optical lattices. The Fermions would be suitable for storage 
of quantum information due to their non-interacting behavior, 
whereas the Bosons could be used to let the systems interact 
and perform operations. 

In turn, the study of such coupled quantum systems exhibit- 
ing collective phenomena with the methods of the theory of 
multi-particle entanglement is an attractive investigation in its 
own right. The aim here is not to use the coupled system to 
prepare strongly entangled systems that form the starting point 
for applications in quantum information processing. Instead, 
the motivation of such investigations is to go beyond conven- 



tional methods to characterize the natural correlations present 
in the distributed system at zero temperature I43ll44ll45ll4fll . 

The quantitative theory of entanglement that abstracts from 
the actual physical realization of a Bose-Fermi Hubbard 
model could provide the tools to understand how global prop- 
erties emerge here from quantum correlations between the el- 
ementary constituents. 



VII. ACKNOWLEDGMENTS 

We thank Kai Bongs and Giovanni Modugno for very use- 
ful discussions. A. A. and J. E. thank the DFG and the ESF for 
financial support. F. 1. thanks the INFM for financial support 
as weU as COSLAB and BEC2000H- ESF programs. 



[1] M. Greiner, O. Mandel, T. E.sslinger, T. W. Hansch, and I. 

Bloch, Nature 415, 39 (2002). 
[2] M. Greiner, I. Bloch, O. Mandel, T. W. Hansch, and T. 

Esslinger, Phys. Rev. Lett. 87, 160405 (2001). 
[3] C. Orzel, A. K. Tuchman, M. L. Fenselau, M. Yasuda, and M. 

A. Kasevich, Science 291, 2386 (2001). 
[4] J. R. Anglin, and W. Ketterle, Nature 416, 211 (2002). 
[5] P. S. Jessen, and I. H. Deutsch, Adv. At. Mol. Opt. Phys. 37, 95 

(1996). 

[6] S. Sachdev, Quantum Phase Transitions (Cambridge University 

Press, Cambridge, 1999). 
[7] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P ZoUer, 

Phys. Rev. Lett. 81, 3108 (1998). 
[8] D. van Oosten, P. van der Straten, and H. T. C. Stoof, Phys. Rev. 

A 63, 053601 (2001). 
[9] J. Ruostekoski, G. V. Dunne, and J. Javanainen, Phys. Rev. Lett. 
88, 180401 (2002). 
[10] W. Hofstetter, J. 1. Cirac, P ZoUer, E. Demler, and M. D. Lukin, 

Phys. Rev. Lett. 89, 220407 (2002). 
[11] B. Paredes and J. 1. Cirac, Phys. Rev. Lett. 90, 150402 (2003). 
[12] A. Recati, P. O. Fedichev, W. Zwerger, and P. ZoUer, Phys. Rev. 

Lett. 90, 020401 (2003). 
[13] H. P Buchler, G. Blatter, and W. Zwerger, Phys. Rev. Lett. 90, 
130401 (2003). 

[14] A. J. Kerman, V. Vuletic, C. Chin, and S. Chu, Phys. Rev. Lett. 
84, 439 (2000). 

[15] P S. Jessen, D. L. Haycock, G. Klose, G. A. Smith, 1. H. 

Deutsch, and G. K. Brennen, Quant. Inf. Comp. 1, 20 (2001). 
[16] L.-M. Duan, E. Demler, and M. D. Lukin, LANL Preprint 

cond-mat/0210564 (2002). 
[17] 1. H. Deutsch, G. K. Brennen, and P. S. Jessen, Fortsc. der 

Physik48 925 (2000). 
[18] D. Jaksch, H.-J. Briegel, J. 1. Cirac, C. W. Gardiner, and P 

ZoUer, Phys. Rev. Lett. 82, 1975 (1999). 
[19] J. J. Garcia-Ripoll and J. 1. Cirac, Phys. Rev. Lett. 90, 127902 

(2003). 

[20] U. Domer, P. Fedichev, D. Jaksch, M. Lewenstein, and P. ZoUer, 

LANL Preprint quant-ph/02 12039 (2002). 
[21] J. Pachos and P L. Knight, LANL Preprint quant-ph/0301084. 

(2003). 

[22] M. P A. Fisher, P B. Weichman, G. Grinstein, and D. S. Fisher, 
Phys. Rev. B 40, 546 (1989). 



[23] For a detailed study of the effects of inhomogeneities in Boson- 

Fermion lattice models cfr: A. Albus, M. Cramer, J. Eisert, and 

F. lUuminati, in preparation. 
[24] G. Modugno, F. Ferlaino, R. Heidemann, G. Roati, and M. In- 

guscio, LANL Preprint cond-mat/0304242 (2003). 
[25] A. P. Albus, S. A. Gardiner, F. Uluminati, and M. Wilkens, Phys. 

Rev. A 65, 053607 (2002). 
[26] L. Viverit and S. Giorgini, Phys. Rev. A 66, 063604 (2002). 
[27] G.-H. Chen and Y. S. Wu, Phys. Rev. A 67, 013606 (2003). 
[28] L. Viverit, C. J. Pethick, and H. Smith, Phys. Rev. A 61, 053605 

(2000). 

[29] J. K. Freericks and H. Monien, Europhys. Lett. 26, 545 (1994). 
[30] K. Sheshadri, H. R. Krishnamurthy, R. Pandit, and T. V. Ramar- 

ishnan, Europhys. Lett. 22, 257 (1993). 
[31] W. Krauth, M. Caffarel, and J.-P Bouchaud, Phys. Rev. B 45, 

3137 (1991). 

[32] S. Kirkpatrick, C. D. Gelatt, and M. P Vecchi, Science 220, 671 
(1983). 

[33] B. W. Wah and T. Wang, Simulated Annealing with Asymp- 
totic Convergence for Nonlinear Constrained Global Optimiza- 
tion, Proc. Principles and Practice of Constraint Programming 
(Springer, Heidelberg, 1999). 

[34] J. A. Nelder and R. Mead, Comp. Journal 7, 308 (1965). 

[35] R. Roth, Phys. Rev. A 66, 013614 (2002). 

[36] M. A. Cazalilla and A. F Ho, LANL Preprint 
cond-mat/0303550 (2003). 

[37] A. P. Albus, F. Uluminati, and M. Wilkens, LANL Preprint 
cond-mat/021 1060 and Phys. Rev. A 67 (2003), in press. 

[38] G. G. Batrouni, V. Rousseau, R. T. Scalettar, M. Rigol, A. Mu- 
ramatsu, P. J. H. Denteneer, and M. Troyer, Phys. Rev. Lett. 89, 
117203 (2002). 

[39] E. H. Lieb and E Y. Wu, Phys. Rev. Lett. 20, 1445 (1968); ibi- 
dem 21, 192(1968). 

[40] P Vignolo, Z. Akdeniz, and M. P Tosi, LANL Preprint 
cond-mat/0304104 (2003). 

[41] IL P Buchler and G. Blatter, LANL Preprint 
L cond- mat/0304534 (2003). 

[42] W. Dur and H.-J. Briegel, Phys. Rev. Lett. 90, 067901 (2003). 

[43] T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 032110 
(2002). 

[44] J. 1. Latorre, E. Rico, and G. Vidal, LANL Preprint 
quant-ph/0304098 (2003). 



11 



[45] K. Audenaert, J. Eisert, M. B. Plenio, and R. F. Werner, Phys. (2003). 

Rev. A 66, 042327 (2002). 
[46] R Giorda and R Zanardi, LANL Preprint |]quant-ph/03041?T] 



